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ABSTRACT 

The Uranian satellite Miranda presents a high inclination (4.338°) and evidences of 
resurfacing. It is accepted since 20 years (e.g. Tittemore and Wisdom 1989, Malhotra 
and Dermott 1990) that this inclination is due to the past trapping into the 3:1 res- 
onance with Umbriel. These last years there is a renewal of interest for the Uranian 
system since the Hubble Space Telescope permitted the detection of an inner system 
of rings and small embedded satellites, their dynamics being of course ruled by the 
main satellites. For this reason, we here propose to revisit the long-term dynamics 
of Miranda, using modern tools like intensive computing facilities and new chaos in- 
dicators (MEGNO and frequency map analysis). As in the previous studies, we find 
the resonance responsible for the inclination of Miranda and the secondary resonances 
associated, likely to have stopped the rise of Miranda's inclination at 4.5°. Moreover, 
we get other trajectories in which this inclination reaches 7°. We also propose an an- 
alytical study of the secondary resonances associated, based on the study by Moons 
and Henrard (1993). 

Key words: celestial mechanics - planets and satellites: dynamical evolution and 
stability - planets and satellites: individual: Miranda 



• 1 INTRODUCTION 

In the 1980s, the Voyager 2 spacecraft gave us a better 
■ know ledge of Uranus and its satellites (see e.g. ISmith et al.l 
' 1 19861 ). It revealed in particular that Miranda and Ariel 
, have been resurfaced. Moreover, orbital models compared 
' to astrometric observations showed that Miranda had a 
significant inclination im, of the order of 4.338° among 
'. GUST8d^ (|Laskaill 19861 ). These facts induced several dynam- 
ical studies in the late 1980s and early 1990s (iDermott et al. 
1988'; 'Malhotra" 1988'; Tittemore & Wisdom 1989; Malhotra 
1990, ; Tittemore & Wisdom 1990; Henrard & Sato 1990) 
showing that the current inclination of Miranda is probably 
due to a former 3:1 resonance with Umbriel. Indeed, tidal 
interactions with Uranus are supposed to push the satel- 
lites outward, meeting orbital resonances. Once the system 
is trapped into this 3:1 resonance, the inclination of Miranda 
is pumped and the amplitude of libration of the resonant ar- 
gument rises as the trajectory meets several secondary res- 
onances. The capture in one of these resonances leads the 
trajectory to the edge of the primary resonance involving 
the exit of this latter. Eventually, the inclination of Miranda 
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ceases to be pumped and the 2 satellites restart to migrate 
independently of each other. This scenario of evolution has 
been intensively studied by the above authors. With this 
work, we confirm the ideas developed in the past and intro- 
duce modern numerical tools to update the problem. 

The reason is that there is a renewal of interest these 
last years for the Uranian system. First, the Hubble Space 
Telescope allowed the discove ry of a whole system of ring s 
and inner satellites (see e.g. IShowalter fc Lissauerl [ioOSi ). 
whose dynamics is of course widely influenced by the main 
satellites. These inner satellites present interesting dynami- 
cal configurations and myster ies, as for instance the poorly 
understood dynamics of Mab (|Kumar et al.|[201ll ) or the in- 
stability of Cupid, Belinda, Cr essida and Desdemona on a 
time scale of 10^ - 10^ years jPrench fc ShowalteJ [ioT^ ) . 
Secondly, a ne w scen ario has recently been proposed by 
iBoue &: Laskai) (|2010l ) to explain the huge obliquity of 
Uranus (« 98°). This scenario involves a former giant satel- 
lite, whose gravitational torque was strong enough to tilt 
the planet (MorbideUi et al."2012'). We can also mention the 
work of iDeicnno ct al. (2011) showing that the presence of 
the current main satellites is consistent with the migration 
of Uranus as predicted by the Nice model. Finally, Uranus 
and its satellites are the ta rget of the proposal of a space 
mission Uranus Pathfinder (|Arridge et all 12012^. For all of 
these reasons, we propose to revisit the dynamical history of 
the main Uranian satellites with powerful numerical means 
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and tools, such as frequen cy ana lysis (Taskar'Tgg^ and the 
MEGNO chaos indicator l|Cincott a & Simo 200Q). 

This article is split into the following sections. First the 
key features of the system are summarised (cfr. Section 
We present the system and the mean motion resonance 3:1 
between Miranda and Umbriel that the system should have 
encountered in its history. The details of the initial condi- 
tions used throughout this study are given in order to group 
the information needed to redo the n umerical analysis. Sec- 
ondly, base d primarily on th e work of iTittemore fc WisdomI 
l|l988l l and iMalhotral l|l990l ). we present the full Hamilto- 
nian (cfr. Section [H]) and its averaged form (cfr. Section Q. 
These two sections remind the description of the system via 
Hamiltonian formalism. The introduction of canonical vari- 
ables as Delaunay variables allows us to develop a perturba- 
tive theory and to obtain an Hamiltonian with two degrees 
of freedom (angle-action) by averaging over the fast angles. 
We also introduce the equations related to the tidal effects 
in the two models. 

We then present the new numerical tools implemented 
to study the system. The numerical integrations of the full 
system on a sufficient large time scale were a major problem 
in the previous studies. We perform numerical integrations 
over 1 Myr on the system of Uranus with its five main satel- 
lites and confirm previous results. We use numerical tools 
like the chaos detector MEGNO to represent the global dy- 
namics of the system (cfr. Section [S}. We also improve the 
resolution and the details of the maps by studying the vari- 
ations of the orbital elements involved in the mean motion 
resonance. These scales are introduced in the section [S] and 
the results on the system are summarised in the section [S] 
These new tools allow us to extend previous studies by the 
introduction of new visualisations of the phase plane of the 
considered problem. We show that the combination of nu- 
merical methods like chaos maps and frequency analysis (cfr. 
Section [S)) allows to detect some particular behaviours of 
the system. We find two regions surrounding the center of 
libration where two secondary resonances are superimposed. 
The choice of the initial condition for a trajectory is then 
primordial for the future of the system: the combination of 
the two secondary resonances tends to increase the chaos 
and to result in another scenario than the exit at 4.5° for 
the inclination of Miranda. This is a part of our conclusions 
and perspectives presented in the last section [T] 



2 KEY FEATURES OF THE SYSTEM 

The system of Uranus has five main satellites, from closest 
to farthest with respect to Uranus; Miranda, Ariel, Um- 
briel, Titania and Oberon. Although these satellites are not 
currently locked into orbital resonances, there are many 
clues indicating probable passages through mean motion 
resonances in the past: we observe important resurfacings 
of Miranda and Ariel and some abnormalities in the current 
orbital elements. In view of these orbital elements of the 
satellites, in particular the high inclination of Miranda (cfr. 
Table [TJ , we analyse a mean motion resonance acting on 
the inclinations. The first one encountered by the system 
in the past is the mean motion resonance Miranda-Umbriel 
3:1 defined by 6 possible resonant arguments: 
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where, in the left column, 9i are the resonant arguments for 
the primary resonances with Ai, the mean longitudes, Qi the 
ascending nodes and, zui the pericenters. The indices M and 
U stand respectively for Miranda and Umbriel. The right 
column is the type of the resonance and corresponds to the 
first non-zero term associated with the cosine of the angle 
6i in the perturbative potential (cfr. Eq[3] in the following 
section). 

These six primary resonances are not well separated 
because the oblateness of Uranus, further indicated by 
the parameter J2, is rather small (cfr. Table [3]) which re- 
sults that the isolated resonance theory is only applica- 
ble for small eccentricities and inc l inations ( Dermott et al.l 
1 19881 : ITittemore fc WisdomI 1 19881 : IMalhotral 1 199d ). When 
the inclination of Miranda increases, this classical theory 
breaks down: it has been shown that the proximity of 
other primary resonances implies the exit of that resonance 
(Tittcmorc & Wisdom 1990) due to a commensurability be- 
tween the libration frequency of the resonant argument and 
the circulation frequency of a close primary resonance, in 
other words, due to the passage in a secondary resonance 
zone (ITittemore fc Wisdom|[l990l : lMalhotralll990l ). 

Our main purpose being the study of the high inclina- 
tion of Miranda, we focus on the primary resonance of type 
/If. The mixed primary resonance ImIu implies also the rise 
of the inclination of Umbriel, which is currently in the same 
order of the inclinations of the three other satellites Ariel, 
Titania and Oberon. For this reason, we do not focus on this 
type of resonance. 

For the numerical analysis of the system, we use as ini- 
tial conditions the values given by the JPlE|. AH the pa- 
rameters used are gathered in Table [T] for the mean orbital 
elements of the satellites, in Table [2] for the physical pa- 
rameters of the satellites and, in Table [3] for the physical 
parameters of Uranus. 

The initial conditions in our numerical experiments are 
fixed to the current values, except for the inclination of the 
satellite Miranda and the ratio of semi-major axes of the two 
satellites Miranda and Umbriel. Considering tidal evolution, 
we modify these elements to reproduce the capture in the 
mean motion resonance 3:1 and to study the consequences 
on the system. Since Miranda is the closest main satellite 
of Uranus, we also use its periods to determine the integra- 
tion steps. All these points are detailed in each numerical 
experiment presented in the following sections. 



3 PI : FULL PROBLEM 

The full problem i.e., the N-body problem with the gravi- 
tational perturbations of each body in the system, the ef- 
fect of the oblateness of the planet and the tidal effect on 
the semi-major axes and eccentricities is described in this 
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Table 1. Mean orbital elements of the five main satellites llLaakar Sz Jacobsonlll98^ : a is the semi- major 
axis, e the eccentricity, ui the pericenter, M the mean anomaly, i the inclination, Q the ascending node, n 
the mean motion. The variables P and Pq stand for the orbital and the node periods respectively. 
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(deg) 


(deg) 
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(deg/day) 


(days) 


(yr) 


Miranda 


129 900 


0.0013 


68.312 


311.330 


4.338 


326.438 


254.6906576 


2.520 


17.727 


Ariel 


190 900 


0.0012 


115.349 
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22.394 
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57.248 
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266 000 


0.0039 
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Titania 


436 300 


0.0011 


284.400 


24.614 


0.079 


99.771 
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13.46 


195.369 


Oberon 


583 500 


0.0014 


104.400 


283.088 


0.068 


279.771 


26.7394888 


1.413 


195.37 



Table 3. Physical parameters and corresponding incertainties of Uranus: QM is given by'jacobson' 1*20071). 
The parameters Re and Rp stand for the equatorial ( Jacobson 2007) and the mean radius (Seidclman c t al.l 
1^07) of the planet re spectively and, J2 and J4 for the spherical harmonics associated with the oblateness 
of the planet ll Jacobson. ,2007, ) . 
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Uranus 5 793 964 ± 6 


26 200 


25 362 


3 341.29 ± 0.72 


-30.44 ± 1.02 



Table 2. Physical parameters and corr esponding incert ainties of 
the five main satellites: QM is given byLj acobson' ('2007') and the 
mean radius of the satellites R by I Thomas 11988, ). 
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R 




(km^/s^) 


(km) 


Miranda 


4.4 ±0.4 


235.8 ± 0.7 


Ariel 


86.4 ±5.0 


578.9 ± 0.6 


Umbriel 


81.5 ± 5.0 


584.7 ±2.8 


Titania 


228.2 ± 5.0 


788.9 ± 1.8 


Oberon 


192.4 ± 7.0 


761.4 ±2.6 



section and will be denoted by PI in the further sections. 
We use a planetocentric reference frame, and consider the 
perturbations of A'^ satellites seen as point masses, and the 
spherical harmonics J2 and J4 of the gravity field of Uranus. 
In this framework, the equations of the problem are: 



Fp 
M ' 



(1) 



where n — {xi,yi,Zi) locates the satellite i, rrn being its 
mass and Fi the force acting on it, Fp and M are respec- 
tively the force acting on Uranus and its mass. We write the 
general equations of motion for the A'^ satellites: 

^ _ g{M + m,) 



connected to Uranus, and P„ the classical Legendre poly- 
nomial. The gravitational potential ([3)l only considers the 
known spherical harmonics J2 and J4 for the Uranian Sys- 
tem (cfr. Table O. 

By introducing Jacobian coordinates, we can write the 
usual Hamiltonian to the first order on satellite masses 
IITittemore fc WisdomlHosil ): 
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where a; is the semi-major axis of the satellite i and TZij 
the disturbing function of the satellite j acting on i. This 
complete Hamiltonian Q considers the mutual gravitational 
interactions between the A'' satellites of the system as well 
as the oblateness of Uranus. 

As dissipation effect, we add the tidal effect on the ec- 
centricities and semi- major axes on the sat ellite i via Kaula 
formulations (see e.g. lYoder fc Peal3ll98ll ): 
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Q being the gravitational constant. Re the equatorial radius 
of the planet, (fn the latitude of the satellite i in a frame 



Rp being the mean radius of Uranus {Rp 7^ -Re), Ri and 
Ui representing respectively the mean radius and the mean 
motion of the satellite i. We observe that these formula- 
tions depend on the Love number ^2 and on the dissipa- 
tion function Q. These are secular equations assuming that 
the satellites are in sy nchronous rotation, a s expected from 
their tidal despinning l|Gladman et al.lll996l ). Moreover, due 
to our poor knowledge of the relevant values, the dissipa- 
tion functions are assumed to be constant with respect to 
the tidal frequencies. 
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4 P2 : AVERAGED PROBLEM 

In this section, we introduce an averaged version of the 
problem PI that will be denoted by P2 in the further 
sections. To reduce time computation and to elaborate an 
analytical perturbative the ory of the problem, following 
iTittemore &: WisdomI l|l988l 'l. we perform an analytical av- 
eraging process on the short period angles Ai. As at the 
lowest order, there is no coupling between the resonances 
in eccentricity and in inclination, the averaged model con- 
siders a cir cular-inclined approximatio n for the inclination 
resonance (|Tittemore fc WisdomI Il989l ) which allows us to 
write the Hamiltonian Q like: 



kep 



Hob + Hre 



(7) 



splitting into the keplerian, the effects of the oblateness of 
the planet parts, the resonant Hamiltonian, and the secular 
one. We obtain finally a Hamiltonian with two degrees of 
freedom ( Jm, Ju, 6m , du), these variables being canonically 
conjugate l|Malhotralll990l ): 
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6i being the resonant angles for the resonance Miranda- 
Umbriel 3:1 described in Section [2] and Ji, the associated 
variables given to th e lowe st order in rui and Si = sin i /; . 
Following iMalhotral l|l988h . we define the parameters de- 
pending on the problem: 
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The expression i^o = corresponds to the exact 3:1 com- 
mensurability between the mean motions of Miranda and 



Umbriel, and the Aui are the corrections of the secular pre- 
cession rates of the nodes Qi on the resonant combination 
of the mean motions of the satellites. The terms b'f\a) are 
the Laplace coefficients, with a = aM/au- Their expressions 
hav e been numerically comput ed from the integral given by 
e.g. iMurrav fc Dermotd (j 19991 ): 

Kic^-ir r^^^'^ ... • (18) 
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We also define the ex pressions of d epending on the inch- 
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where we neglect the indirect perturbation of order niM/M. 

To follow a trajectory through the resonance 3:1 in in- 
clination with the tidal effects, we also consider the Kaula's 
formulations © where we assume the circular approxima- 
tion. We obtain a tidal perturbation on the semi-major axis 
of both satellites Miranda and Umbriel written as: 

dt Q^afM ' 

where the dissipation factor of the satellites (k^/Qi) is ne- 
glected in the case of the averaged problem. Because of the 
circular approximation, we have a tidal perturbation only 
on the semi-major axes of the satellites. 



5 NUMERICAL METHODS 

A major obstacle in previous studies was due to the numeri- 
cal integrations of the full equations of PI on a long enough 
time scale. The rate of the tidal evolution of the system 
being rather slow, the CPU time needed is quite high. 

Using new computing tool^, we have integrated nu- 
merically the equations of motion of the N-body prob- 
lem in cartesian coordinates (cfr. Eq[2]) with the Adams- 
Bas hforth-Moulton 10* ^ order predictor-corrector integra- 
tor ()Hairer et al.l l200ll ). and the outputs have been suc- 
cessfuUy compared to t he well-known software SWIFT 
(jLevison fc Duncanlll994 ) on the Uranian system. 

^ The computations were performed on an HPC cluster which 
powers the Interuniversity Scientific Computing Facility cen- 
ter located at the University of Namur (Belgium). iSCF - 
|http://www.iscf.be| 
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In this section, we introduce several methods used to 
obtain our numerical results. We first present the chaos in- 
dicator call ed Mean Exponential G rowth of Nearby Orbits 
(MEGNO) (|Cincotta fc Simdl2000l ). We perform a set of nu- 
merical integrations on the problems PI and P2 to represent 
the phase plane of the problem with a chaos map. 

Secondly, we study the variations of the orbital elements 
of the satellite Miranda. According to the type of variations, 
different details of the eye of the resonance 3:1 Miranda- 
Umbriel, which are not visible with the chaos indicator, ap- 
pear. More information is following in the subsection 15.21 



5.1 Chaos indicator 

To study the stability and the structure of the phase space of 
our system, we implement the chaos indicator called Mean 
Exponential Gr owth of Ne arby Orbits (MEGNO) elabo- 
rated by iCincotta fc Sim^l (2000) and used by different 
authors (i.e. iGozdziewski et all i 200ll : ICincotta etal] l2003l : 
IValk et al.|[2009l : lDelsatell201ll ). The indicator MEGNO is 
defined by: 



Y{t) 



m 



(28) 



where 5 — \ \X\\, X being the tangent vector {Sp,Sq) of the 
vector angle action (p, q) and 5, the time derivative of 5. 
Given this last equation (|28|l . the running time-average of 
the MEGNO, Y{t), is defined by the foUowing integral 



Y{t)^\j\{t') dt', 



(29) 



where, as explained by ICincotta et al.l l|2003h . Y{t) 
quickly the behaviour of the system, i.e. 



lim Y{t) = 0: stable and periodic trajectory , 
lim Y{t) 



2: stable and quasi-periodic trajectory 



hm^ 

t— >oo t 



> 0: chaotic trajectory . 



IGozdziewski et all l|200lD propose an efficient way to 
compute the MEGNO, by adding two differential equations 
to the system: 



dy 
dt 
dw 
HI 



5-5 
5- 5 
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(30) 
(31) 



and by computing, at each step of the integration, the 
MEGNO and its averaged by: 



Y(t) = 
Y{t) = 



t 

t 



(32) 
(33) 



In a first step, integrating the equations (|30|) and (1311) . 
we compute the indicator using the definitions (|32|l and (|33p 
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Resonant Argument S ^=X^-5X^+2Q^ (deg) 

Figure 1. Chaos map in the plane semi- major axis ajv/ vs. res- 
onant argument 9m resulting from the numerical integrations of 
the 3 body problem Uranus-Miranda-Umbriel with the Adam- 
Bashforth-Moulton integrator over 1 500 years. The integration 
step is setted to 1/80 day. The initial conditions are fixed to the 
current ones (cfr. Tables [T]|3] and [2)l except for the mean anomaly 
Mm, the semi-major axis aM and the inclination im of the satel- 
lite Miranda. The two first variables are set respectively between 
[0° - 360° [ and [127845 km - 127895 km]. The initial inclination 
is fixed to 4.8°. 



on the problem PI with Uranus and its five main satellites 
perturbing each other. In a second step, again on the prob- 
lem PI, we compute the indicator by considering a 3 body 
problem consisting of Uranus, Miranda and Umbriel, these 
two satellites being involved in the mean motion resonance 
studied. In both cases, the first deviation vector is chosen 
randomly. We notice that the main points of the global dy- 
namics are similar in the two simulations and conclude that 
considering a 3 body problem is sufficient to represent the 
evolution of the system during the passage in the 3:1 reso- 
nance. 

To obtain the chaos map (cfr. Figure [1} in the plane 
semi-major axis aj\/ vs. resonant argument 9m, we perform 
10^ numerical integrations over 1 500 years, each one asso- 
ciated with an initial condition of the phase plane. We fix 
the integration step at around 1 /SO*'' of the smallest period 
of the five main satellites (cfr. Table [l]). For each numeri- 
cal simulation, we compute the solutions of the equations of 
motion of the 3 body problem PI, without adding the tidal 
perturbation, and the averaged MEGNO value (|33|l to have 
the third dimension for the colour scale. The tidal effect is 
not considered because instead of following one trajectory 
in time, we select a region of initial conditions which covers 
the entire dynamics (um G [127 845 km — 127 895 km] and 
Mm G [0° —360° [) . We do not need the effect of tides to push 
the pair of satellites in the resonance. This method allows 
us to have a global visualisation of the eye of the resonance. 

In Figure [T] we obtain the phase plane given by the av- 
eraged MEGNO value. The indicator identifies the chaotic 
and the stable structures of the phase plane (resp. dark and 
light colours): it detects the external separatrix but also a 
smaller one defining the boundary of secondary resonance 
zones. Indeed, we note first the wide separatrix between the 
libration and the circulation zones. The stable zone in the 
large separatrix (top right in the figure) is the following pri- 
mary resonance in inclination considering the nodes of MI- 
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randa Qm and Umbriel ilu, represented by the resonant 
argument 82. The other structures appearing in the circu- 
lation zones are the primary resonances 3:1 in eccentricities 
9i, with i = 4, 5,6. The visualisation of all these structures 
in a small range of semi-major axis confirms the well-known 
narrowness between the resonances of the system due to the 
small oblateness of the planet. Secondly, we distinguish an- 
other separatrix in the libration zone suggesting the presence 
of secondary resonances zones (in the middle of the figure) 
that will be examined more closely in the following sections. 



5.2 Orbital elements variations 

The second method in our numerical analysis considers the 
variation in orbital elements of the satellite Miranda. When 
the bodies pass through a resonance, they cross the sepa- 
ratrix between the circulation and libration zones inducing 
changes in the nature of the orbit: we observe modifications 
in eccentricities or inclinations of the satellites involved in 
the resonance. In contrast, in the center of the resonance, 
the variations are negligible as the orbital elements are con- 
stant. Indeed, the closer is the separatrix, the stronger is 
the change. We base our new colour scales, named orbital 
element variations, on these previous facts and define the 
following colour scales for the map: 

(i) the variation of the semi-major axis of Miranda, 
Sa (km) , 

(ii) the variation of the inclination of Miranda, 5i (deg) , 

(iii) the variation of the eccentricity of Miranda, 8e ■ 

The study of the variations of the orbital elements allows 
the visualisation of structures of the space phase invisible 
with the chaos indicator MEGNO: these structures are 
not chaotic and/or are too barely perceptible because the 
variations are very small. We also have the association of the 
structures with the orbital element involved. To illustrate 
this, we consider the full problem PI. and represent the 
phase plane of the 3:1 resonance between Miranda and 
Umbriel (cfr. Figure [2| obtained by the same integration 
process than the Figure[TJ but with the inclination (cfr. Fig- 
ure[2](a)) and the eccentricity (cfr. Figure[5] (b)) variations 
for the colour scale. What we observe is the following : we 
have the eye of the resonance with the inclination variation 
for the colour scale whereas it is almost erased in the 
eccentricity one. Moreover, with the eccentricity variation 
for the colour scale, we have another structure that does 
not appear in the inclination variation map: this is another 
primary resonance involving the pericenter of Miranda -cum . 
The semi-major axis variation map brings together the 
structures of the inclination and the eccentricity variations 
maps giving a global visualisation of the space plane (cfr. 



1 .2790 



0.15 



Figure [6} . 



6 RESULTS AND COMMENTS 

We apply the numerical methods presented in the section [S] 
to the 3 body problem Uranus-Miranda-Umbriel to study 
the 3:1 mean motion commeiisurability between Miranda 



0.05 




100 200 300 

Resonant Argument %^=X^-?,X^+2Q.^ (deg) 

Figure 2. Eyes of the resonance from the 3 body problem 
Uranus, Miranda, Umbriel. The integrator, the integration step, 
the model, the initial conditions are the same as Figure[T] The ini- 
tial inclination of Miranda im is fixed to 4.338° . The colour scales 
consider the variations in the inclination of Miranda (deg) (a) 
or the variation in eccentricity of Miranda ejv/ (b). We see the 
structures related to the resonance in inclination involving the 
node of Miranda (a) and, the structures related to the resonance 
in eccentricity involving the pericenter of Miranda ro^/ (b). 



and Umbriel. We consider both models PI and P2, the lat- 
ter allowing us to speed up the CPU time while maintaining 
the global evolution of the system. 

We have two different types of results. The first one 
groups the trajectories which evolve in the time with a tidal 
dissipation on a large time scale, typically 1 Myr. The usual 
problem in this case is the dissipation function Q which 
value is not known but recent studies suggest s maller values 
for p lanet dissipation factor than ever (see e.g. iLainev et al.l 
12012! for Saturn) . We implement a set of tests with different 
values for k2/Q and we note similar behaviours on different 
time scales. We choose the value of 5.2 10~^ for the planet 
because it produces the rise of inclination of Miranda signif- 
icantly higher than the current one, allowing us to study the 
overall process in a moderate CPU time. The strengthening 
of the parameters k2/Q allows us to speed up the integra- 
tions provided that the trajectories remain adiabatic. 

Secondly, we represent the global evolution via maps de- 
pending on a colour scale MEGNO to represent the separa- 
trix or depending on orbital elements variations, principally 
the inclination one, to study with more details the resonance. 
In this C£LS6j £LS 3, reminder, we represent the phase plane in a 
domain semi-major axis vs. resonant argument determined 
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by a set of numerical integrations over 1 500 years without 
taking into account the tidal effect. 



6.1 Type of results 1 : the tracking of a trajectory 

Let us consider the tracking of a trajectory in its capture 
in the resonance 3:1 between Miranda and Umbriel. The 
equations of motion derive from the model PI (cfr. Section 
[S]). At the begining of the simulation, the initial orbital el- 
ements are set to the current ones (cfr. Table [T]) except for 
the semi-major axes of Miranda and Umbriel and the in- 
clination of Miranda which are fixed at lower values. We 
add the tidal effect on the semi-major axes and eccentrici- 
ties of Miranda and Umbriel following the Kaula's formula- 
tions ((BJ. The ratios k2/Qi are fixed to 5.2 10~^ for Uranus, 
10^* for the satellites. We can also assume a dissipation fac- 
tor Q higher for Mirand a than Umbr i el as Miranda shows 
an important resurfacing (jSmith et al.l (|l986l )) suggesting an 
important past thermal history but this kind of aspect will 
be study soon in a following paper. 

The result of the numerical integration gives us the evo- 
lution of the resonant argument 6m during the capture in 
the resonance. In the Figure [31 the tidal evolution of the 
resonant argument is plotted vs. the semi-major axis ratio 
a,M /au vs. the rise of inclination of Miranda ja/. We observe 
first a circulation of the angle ^a/. The capture in the res- 
onance Im occurs when the commensurability 1/3 between 
the semi-major axis is approaching, corresponding to a ratio 



a A/ 
au 



3/2 



= (0.4807)^''^ 
= 0.3333 . 



Then we enter in a libration zone and the inclination starts 
to increase. The trajectory is later captured in a secondary 
resonance zone which implies the exit of the primary reso- 
nance zone. Let us note the changes of regime at the same 
critical points than .Malhotra (jl990i ) (ia/ = 2.1°, im = 2.8° 
and lAf = 4.1°). 



6.2 Type of results 2 : Global visualisations via 
phase planes 

6.2.1 Averaged problem P2 vs full problem PI 

Let us consider the averaged system consisting of a 3 body 
problem Uranus, Miranda and Umbriel in a circular-inclined 
approximation described by the model P2 (Cfr. Section U). 
To compare the quality of the averaged problem P2 vs. the 
full one PI we implement the same process explained in the 
subsection l5.2l on the two problems. The integration step for 
the problem P2 is fixed to 17/300 year, i.e. 1/300"^ of the 
smallest nodal period. In our case it is the period of 17 years 
related to the node of Miranda (cfr. Table [T|. The Figure 
[4] compares two eyes of the resonance obtained by numeri- 
cal integrations of the equations of motion of the complete 
model PI for the first one (a), and the equations of motion 
of the averaged model P2 for the second one (b) with the 
inclination variation in degrees for the colour scale. 

We observe that as a whole, the two eyes are similar: we 
distinguish the libration and the circulation zones, and the 
stable zone for the mixed resonance 02 also appears in the 



two models. We observe also the two zones of secondary res- 
onances confirming the possibility of their studies with the 
averaged problem P2. Obviously the details and the preci- 
sions of the eye in the complete model PI are not stored 
in the averaged problem P2 but the global dynamics is well 
represented with the approximation. We note that the colour 
scale is identical. It is due to the fact that the dynamics of 
the inclination is ruled by the resonant (or quasi-resonant 
out of the separatrix) argument, present in the two systems 
(PI and P2). To check our assumption we use the NAFF 
algorithm as Numerical Analysis of the Fundamental Fre- 
quencies based on Laskar's original idea (s e e for instance 
iLaskaj (119931 ) for the method, and iLaskaij l|2005h for the 
convergence proofs). It aims at identifying the coefficients 
Ok and uj^: of a complex signal f{t) obtained numerically 
over a finite time span [— T; T] and verifying 



fit) ~ ^fflfcexp(itjfet). 



(34) 



where LUk are real frequencies and complex coefficients. 
Using this tool, we check that the variations of the inclina- 
tion have the same period as this resonant angle in the two 
problems PI and P2. 

A last thing in the comparison of the two problems is 
the following: looking at the two Figures [4] (a) and (b), we 
notice a shift of the eye of the resonance in the two prob- 
lems. Indeed the two ranges of semi-major axis for Miranda 
are not the same in the two problems. This shift is due to 
the definitions of the variables Jm (|lip and Ju (|12p which 
imply that the semi-major axis of Miranda in the averaged 
problem differs by a quantity of order (siniAi/2)^ from the 
complete one. This can be important as the incl i nation in- 
creases. Based on the second chapter in lMalhotral l|l988h . we 
show that: 



3 4) 



+ Jaij a 



+ 3s 



(35) 



where om is the semi-major axis of Miranda in the averaged 
problem and Sj = sinij/2, j being M for Miranda and, U for 
Umbriel. The Figure [5] shows the shift of the eye of the res- 
onance with the increasing of the inclination of Miranda ia/ 
which is significant when the inclination is high. 

6.2.2 Global evolution of the system 

The representations in three dimensions of the phase plane 
(semi-major axis vs. resonant argument with a colour scale) 
is not sufficient to represent the entire dynamics as it evolves 
with the rise of inclination of Miranda. We know that the 
chaotic ity in the system increases w ith the inclination of Mi- 
randa (|Tittemore fc WisdomllTgsi ). As we want to present 
the whole dynamics, we propose to follow the evolution of 
the system with a succesive set of maps in the semi-major 
axis variations for the colour scale resulting from numerical 
integrations of the problem PI. 

The Figure [5] shows a set of six phase planes semi-major 
axis om vs. resonant argument 9i with the variation in semi- 
major axis in colour scale. The initial inclination of Miranda 
is set to 1°, 2.1°, 2.8°, 4.1°, 4.338° and, 4.8° respectively. 
When the inclination increases, the separatrix broadens to 
become a layer of chaotic motion, in particular when two 
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Figure 3. Evolution of the resonant argument 9m vs. the semi-major axis ratio aM/o,u vs. inclination of Miranda im- The initial 
conditions are the current values (cfr. Tables [T1 [3l and [2)l except for the semi-major axis of Umbriel, ajj = 265200 km, the semi-major 
axis of Miranda, ajvf = 127400 km which evolves writh the tidal effect on the semi-major axes and on the eccentricities. The inclination of 
Miranda is fixed to 0.001°. The integration step is fixed to 1/60 day. The ratio k2/Q is fixed to 5.2 10~^ for Uranus, 10"'' for the satellites. 
We note the capture in the primary resonance when the commensurability 1/3 is approached and the consequent rise of inclination of 
Miranda. We also observe the disruption of the primary resonance following the capture of a secondary resonance. 



© 2002 RAS, MNRAS 000, [THH] 



Miranda's dynamical history 9 



(b) iM=2.r 
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(f) iM=4.8'' 
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Resonant argument 9^ =X.|^-3X.u+2£2(^ (deg) 



Figure 6. Phase planes of the resonance from the 3 body problem Uranus, Miranda, Umbriel. The integrator, the integration step, the 
model, the initial conditions are the same as Figure [T] The initial inclinations of Miranda im are fixed to 1°, 2.1°, 2.8°, 4.1°, 4.338° 
and, 4.8° respectively in the Figures (a), (b), (c), (d), (e) and, (f). The third dimension considers the variations in the semi-major axis 
of Miranda aj(/(km). As the inclination of Miranda increases, the separatrix widens approaching the next primary resonances involving 
chaos by overlap. The zones of secondary resonances seem to evolve with the rise of the inclination. 
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Figure 4. Phase spaces from the 3 body problem Uranus, Mi- 
randa, Umbriel for the full problem PI (a) and the averaged 
problem P2. The integration step, the model and the initial con- 
ditions arc the same as Figure [T] in the case of the Figure (a). The 
initial conditions for the averaged problem P2 reflects the shift 
in semi-major axis aM which is set to [128 268 km — 128 318 km] 
and the integration step is fixed to 17/300 day. In both figures, 
the integrator is the same as Figure [T] and the initial inclination 
of Miranda is 4.338°. 



resonances meet. This increase of chaos is due to the overlap 
of two close separatrices, conse quence of the close n ess o f 
the resonances. The Figure 9 in iMoons fc HenrardI (|l993l ') 
shows 'the landscape' of the problem with the location of the 
separatrix of the primary resonance and the centers of the 
secondary resonances in a particular plane. It also shows the 
location of the chaotic layers around the separatrix. We see 
that these layers become larger as the inclination of Miranda 
is high and we observe exactly the same feature in our maps. 

The different zones of the phase planes evolve too and, 
in particular, the center of libration presents different struc- 
tures moving with time: we distinguish some zones of sec- 
ondary resonances which appear and move with the rise 
of inclination. Some of t hem have already been detected 
by iTitteinore fc WisdomI lll990l) and extensive ly studied by 
iMalhotral (|l99d ') and I Moons fc HenrardI l|l993t) . As we were 
intrigued by the different zones in the libration center and 
their evolution with the rise of inclination of Miranda, we 
tried to make a zoom of this center. The Figure [7] shows an 
enlargement of the center of the eye when the initial inclina- 
tion of Miranda is set to 4.338°. We clearly see three zones: 
the center of libration and two other zones of secondary res- 
onances surrounding the center. 
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Figure 5. Averaged semi-major axis of Miranda a^j vs. inclina- 
tion of Miranda ijv/ (deg). The initial condition for the semi-major 
axis of Miranda aM is set in the equation l|35|l to the center of 
libration in the complete model i.e. aM = 127 870 km. The val- 
ues of the semi-major axis of Umbriel ajj and for the inclination 
ifj are the current ones. We see that the shift of the center of 
libration between the two models increases with the inclination 
as a power law. 



6.2.3 The secondary resonance zones 

iDermott et all (j 19881 ) and iTittemore fc WisdomI (|l988h 
show that the role of secondary resonances in the resonance 
3:1 is crucial: a capture into a secondary resonance can ex- 
plain the escape of the primary due to the chaotic layers 
present near the separatrix. To well understand the dynam- 
ics, we identify these secondary resonances via frequency 
analysis and, in particular, the two zones represented in the 
Figure [T] We select an initial condition in one of these zone 
(symbol 'X' in the FigureO and use the resulting trajectory 
for the frequency analysis. We extract the phase of the libra- 
tion in the primary resonance, and combine it with the other 
arguments to reconstruct the argument of the secondary res- 
onance. We write the libration argument Oi : 



cos { 5i t + I 



(36) 



where Ai and 4>i are respectively called the amplitudes and 
the phases of the sinusoidal function. The parameters 5i are 
the libration arguments of the function. We extract then the 
principal eigenmode of the libration A, which is associated 
with the bigger amplitude A, and combine it with the prin- 
cipal eigenmode of circulation of the arguments 62 and 63. 
Let us call them k and x respectively. We obtain the result 
in Figure[7]where we have librations for two different combi- 
nations in a same time: a combination between the libration 
of the resonant argument 9i 

• and the circulation of the argument 62 , 

• and the circulation of the argument ^3 , 

indicating two secondary resonances of type 2/1 and 4/1 
respectively. 

Although they are very similar, they are not the same 
and they present an interesting phenomenon as they ap- 
pear in the same time. We have different elements that indi- 
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Figure 7. Zoom of the center of libration from the 3 body problem Uranus, Miranda, Umbriel (Top). The integrator, the integration step, 
the model, the initial conditions are the same as Figure[T]except for the semi-major axis of Miranda aM S [127860 km — 127880 km]. The 
initial inclination of Miranda im is fixed to 4.338°. The symbol X in the map represents the initial condition for the trajectory analysed 
by frequency analysis (Below). Two different combinations of frequencies seem to librate. We have the first one with the circulation 
frequency of 62 (a) and the second one with the circulation argument 83 (b). 
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cate the difference between these two secondary resonances. 
First, if we change very slightly the initial conditions in the 
dissipative simulation, we can have one combination in a 
libration and the other in a circulation regime or the two 
combinations in a libration regime (cfr. Fig. [8]). 

More rigorously, we separate the two secondary 
resonances by an analytical approach using the simple 
model o n a restricted three body p roblem described by the 
authors iMoons fc HenrardI ()l993l l. The difficulty here is 
to have two comparable frequencies, both in a circulation 
regime. It is usual to use the angle-action variable to 
transform the libration frequency of the primary resonance 
9i in a circulation frequency to be compared with the 
circulation frequen c ies of the other arguments. Following 
[m oons fc HenrardI (| 19931 ). we introduce the canonical 
variables: 



2cr = - Am + 3Ai7 + 2qM S = 
2v ^ \m - 3Xu - 2qu N 



2L, 



2L* + 1 



where 2a — —6\ and 2u — 63. The variables used in these 
definitions are the usual modified Delaunay's elements in 
case of the restricted 3 body problem : 



Li (1 — cosii) 



\i — mean longitude 
Qi — -longitude of the node 

where i stands for M or U for Miranda or Umbriel respec- 
tively. The variable L* is the value of L at the 'exact reso- 
nance'. The Hamiltonian in these canonical variables is de- 
fined by 

H = Ho + Hi , 

= C{N - Sf + AN + 2DScos2a , 

+ diiuV2S cos{a — v) + d2iuV2S cos{a + u) 

+ daiff cos{2a) , (37) 

where the dj a r e co nstant coefficients defined by 
iMoons fc HenrardI (| 19931 ). Introducing the couple angle- 
action (^l/, J) we show that at the equilibrium (a = n/2, 



* = (1 ^ iV)^''^ 



(38) 



where C, D are constant defined in iMoons fc HenrardI 
(|l993l ). We compare this circulation frequency to the cir- 
culation frequency of 



and 
d_ 

dt 



Ms 
dt 

2A + 4C{N - S) 
2A-4D , 



d92 
dt 

oq 

A + 4C(Af - S) 
A-AD . 



(39) 



(40) 



The expressions of t he co nstants are given by the au- 
thors [MoonT^^^enrar^ (|l993l ) and note that the dominant 
term is the constant A. We show that the two different fre- 
quencies are different and can thus quantitatively distinguish 
the two secondary resonances. 



7 CONCLUSIONS AND PERSPECTIVES 

In this work, we focus on the mean motion resonance 3:1 
between Miranda and Umbriel and try to explain the high 
inclination of Miranda. This problem was studied by numer- 
ous authors twenty years ago but the update of some results 
with new numerical tools gives new views of it. We retrieve 
the main results than these authors and improve the under- 
standing of the problem with new powerful methods. 

The chaos detector MEGNO has never been applied in 
the case of the main satellites of Uranus. We show that the 
combination of the chaos detector and the orbital element 
maps brings a new visualisation of the phase planes of the 
problem. The use of maps and frequency analysis on particu- 
lar trajectories allows the detection of unusual zones in these 
phase planes where two secondary resonances are combined 
and dominate the future evolution. 

We show that sometimes, a trajectory could be cap- 
tured in two secondary resonances in a same time: a 2/1 
secondary resonance with the circulation argument of 62 and 
a 4/1 secondary resonance with the circulation argument of 
63, the two being very close. Following the toy model of 
[m oons fc HenrardI l| 19931 ). we derive an analytical expres- 
sion for a circulation frequency ^ for the primary resonant 
argument in libration needed to compare with the other ar- 
guments in circulation. We find two different commensura- 
bilities and can therefore distinguish the two secondary reso- 
nances. The mix of these two secondary resonances increases 
the chaos and leads to another scenario for the evolution of 
the system with an exit of the primary resonance higher 
than 4.5°. 

The dynamical aspects of the Uranian system are nu- 
merous and full of interest. In the future, we will study a 
combination of this dynamical model with an intern evo- 
lution of the satellites. This combination is interesting by 
many points but the main one according to us is the un- 
derstanding of a dynamical abnormality (high inclination of 
Miranda) and a geological anomaly (differentiation of Mi- 
randa) through a resonance phenomenon. 
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as Figure [3l except for mean longitude Aj(/ of Miranda which is increased by a multiple of 9. On the left, we have the evolution of the 
inclination of Miranda im vs. time t associated on the right with the two combinations 2A — k and 4A — x in degrees vs. time t. In the 
first case, we observe the rise of inclination over 4.5° (a) and the libration of the two combinations in a same time (b). In the second 
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